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SUMMARY 

We have developed a theoretical model for the behavior of 
sunspots as individuals. Two models, called the current sheath 
model and the snowplow model, known in the theory of the pinch 
effect are worked out for the present problem of the sunspot 
expansion. Using the observational fact that the magnetic field 
of the sunspots grows to about 3,000 gauss in about ten days, 
numerical calculations based upon the current sheath model show 
that the sunspot area grows with the magnetic field and begins 
to pulsate when the magnetic field stops growing. The amplitude 
and frequency of oscillation depend upon the mass in the current 
sheath and the maximum magnetic field. The model suggests 
that the bright ring about the sunspots may exist at chromo- 
spheric height where the density of the material is about 
10 9 particles/cm 3 . 
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A THEORETICAL MODEL FOR SUNSPOT COOLNESS 


by 

R. K. Jaggi 

Goddard Space Flight Center 


INTRODUCTION 

Among the unsolved problems in solar physics reviewed by Goldberg and Dyer (Reference 1) is 
the difficulty in obtaining a theoretical interpretation of the mechanism of sunspot cooling. They 
proposed two possible explanations: The gas pressure in the sunspot is low because part of the total 
pressure is exerted by magnetic forces; in that case the configuration must approach an equilibrium 
in which the sum of the kinetic and magnetic pressures equals the kinetic pressure outside. The 
cooling could also be due to the forced expansion of a rising gas column. This second possibility is 
discarded because there is no observational evidence that the gas column is really rising. Many 
other models used to explain the coolness are based on the presence of a strong magnetic field (Ref- 
erence 2, p. 172), Biermann (Reference 3) assumed that the low temperature is probably maintained 
by the strong magnetic fields inhibiting convective transport of energy, thus causing steep tempera- 
ture gradients in the outer parts of the spots. 

At this point it would be instructive to estimate the time in which the sunspot region, radiating 
like a blackbody, would fall from 6500 to 4500°K, the temperatures of the photosphere and a typical 
sunspot respectively. If convection is completely stopped, the time of this temperature fall is deter- 
mined by 


L -gif (nkT) - aT 4 , 


where a - 5.6 x 10" 5 erg/cm 2 -sec-deg 4 , and L , the total depth of the cool region, is given by 


L 


nk( 6500 - 4500) 

- 


ct( 5000 ) 4 , 


with k the Boltzmann constant and a the radiation constant. For L = 10 8 cm and n = 10 17 , the 
above equation gives t £ 7.8 sec, which shows that the region will cool down instantaneously. Some 
work has been done on the inhibition of convection in the presence of a magnetic field (References 4 
and 5). The difficulty with BiermamTs explanation is twofold: (1) it fails to explain the growth of 

the sunspot area with time; and (2) the cooling appears one or two days after the observation of the 
magnetic field. Also, any explanation of the sunspot coolness must at the same time explain the 
Ever shed effect as well as the bright ring around the sunspots (as known by observation). 
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De Jager (Reference 2) has mentioned that the ions with energies greater than 3ev are reflected 
by the sunspot's magnetic field; but this is only effective in the upper, less dense part of the solar 
atmosphere and therefore does not explain lower temperature at the photospheric level. The parti- 
cles reflected by the strong magnetic field of the sunspot will leave the region only if the mean free 
path is large enough to permit their escape. From Table 1 (in the last section of this paper), it is 
clear that the mean free path is large only in the upper chromosphere. 


OBSERVATIONAL DATA 


Before considering the theoretical model, we shall present some observational data. Most of 
this material was given in two articles by De Jager (Reference 2 and 6). 

The development of sunspots forms part of the development of a center of activity (CA). 

De Jager (Reference 6) has divided the development of a CA into four parts: (1) pre-spot phase; 

(2) spot phase; (3) bipolar magnetic post-spot phase; and (4) unipolar magnetic phase. Of these, the 
second is of most interest in the present problem. A typical sunspot consists of an umbra and a 
penumbra, where the latter shows a radial filamentary structure, being slightly brighter than 
the former. The mean umbral area of large sunspots is of the order of 5 x 1CT 4 of the sun's visible 
hemisphere and the radius of the umbra is of the order of 2.2 x 10 9 cm. 

In the pre-spot phase the development of a bipolar magnetic field region as well as faculae 
are observable without the spots. The limits of the faculae coincide more or less with the sunspot 
region. The spot appears one or two days after the appearance of the magnetic field. The spot and 
the magnetic field region increase in area simultaneously. 

A fifteen day time sequence of the development of a typical sunspot group is given by De Jager 
(Reference 2, p. 167). On the first day a single spot or a spot group without penumbra or without 
bipolar structure is observed on the sun. On the second day the bipolar group and penumbras of 
some of the bipolar spots appear. From five to ten days great bipolar spots as well as many small ones 
appear. From ten to fifteen days great bipolar spot groups remain without their smaller companions. 
During this time unipolar spots may also be observable. 

Kiepenheuer (Reference 7) quotes Thiessen (Reference 8) as giving the following radial dis- 
tribution of the vertical magnetic field in a sunspot: 


H,(r) 



( 1 ) 


where b is the radius of the outer edge of the penumbra, the subscript z denotes the z component 
of the field and the subscript m the field at the center of the spot, and r the radial distance meas- 
ured from the axis of the spot. However, the magnetic field distribution as determined by Mattig 
(Reference 9) is of the form: 


«,(r) 
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Figure 1 is an approximate representa- 
tion of this field (taken from Reference 2, 
p. 156). At the center of the spot the 
magnetic field H m is of the order of a few 
thousand gauss. Assuming H m to be 10 3 
gauss we find that the magnetic pressure 
h 2 / 8 t 7 is approximately 4 x 10 ’ and is 
much higher than the kinetic pressure 
of the solar photospheric plasma nkT £ 
3 x 10 3 . Thus the sunspot is under the 
action of strong magnetic forces. Equating 
H 2 / 8 v with the kinetic pressure 3 X 10 3 , 
we obtain H a 260 gauss. Therefore 
whenever the magnetic field is compara- 
ble to or higher than this value the sun- 
spot will be under the action of strong 
magnetic forces. 



Figure 1— Variations of magnetic field strength over a 
sunspot as function of radius from the center of the spot. 
The distance to the spot center is expressed in terms of 
penumbral radius. 


The time variation of the magnetic 
field observed over the sunspots is pre- 
sented in Figure 2, which is taken from 
Cowling (Reference 10). Analytically the 
curve may be approximated by 

H m - 3000^1 - e" 2 ^1 - e 5 - 5 (2) 



where t is measured in days. Figure 2-Variation of the field strength versus time in 

days for a 55-day spot. 

We shall now discuss the adiabatic 

expansion of sunspots for the two models known as the current sheath and the snow plow models, in 
order to determine the rate of cooling in both instances. 


THE CURRENT SHEATH MODEL 


The basic equations of magnetohydrodynamics are 


dv 1 . 

Nm a "dt r T J xK ' K rad P 


3N 


TVf “ -div(Nv) 


Curl B pr j 


(3) 


(4) 
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div B 


c)B 

Curl ( v x B) , 


div E 477 1 , 


where the symbols have their usual meaning, c is the charge density, and m. the mass of an ion. 
From Equations 1 and 4, we obtain, in cylindrical coordinates 


therefore the force, 




Figure 3— Modified form of Figure 1 used in calculating 
the expansion of sunspots. 


is directed radially outward force. Consequently 
the cylinder must expand in the radial direction. 

To simplify the model we modify the 
distribution of the magnetic field given by 
Equation 1 in Figure 1 to that given in Fig- 
ure 3. A portion of the cylinder is assumed to 
have a uniform magnetic field while the gradients 
of the magnetic field are assumed to be in the 
rest of the cylinder. The region occupied by this 
current is called "the current sheath." Equation 
(3) can be rewritten as 


Nm. 


dv 

dt 


877 


grad B 2 + 4 ^ (B • V ) B - grad p 


In the present case the term (b • v) B = 0 because it is assumed that b is in the z-direction but there 
are no gradients in the z-direction. Integrating over the thickness of the current sheath, multiplying 
both sides by 2?7r , assuming that the acceleration over the cross-section of the current sheath and 
the density in the current sheath are uniform, we obtain 


M d 2 r 
a dt 2 


/V 

27Tr l 877 + Pi 



(5) 


where p 0 is the kinetic pressure due to the gases outside the cylinder, p. the average kinetic pres- 
sure due to gases within the cylinder, and B i the average magnetic field within the cylinder, and M 
the total mass contained in the cylinder, a has been introduced to insure that only a fraction of the 
total mass of the cylinder is contained in the current sheath. We shall assume that p 0 remains 
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constant and that p. = Pfl at t = 0. The changes in the kinetic pressure satisfy the ordinary adia- 
batic law 


— “ ; Constant . (6) 

N 5/3 

If B 0 is the initial magnetic field and r 0 the initial radius of the cylinder, the constancy of the mag- 
netic flux in the cylinder gives 


77 r 


o 2B o 


B. , 


or 



Then, since N is proportional to 1/r 2 , Equation 6 yields 

- - M 10 ' 

P, • Po^r j 


(7) 


( 8 ) 


Substituting Equations 7 and 8 into Equation 5 gives 


M d 2 r 

2rui 2 




Now by writing r/r 0 = x , B$/8mp 0 


1 

and (m/2vnp 0 ) 2 t' - t the above equation becomes 


1 d 2 x _ d 1 (9) 

* dt' 2 x< + if " ‘ 

To discuss the solution of Equation 9, the dependence of B on time must be known. According to 
observations, Equation (2), and Figure 2, the magnetic field rises in a matter of 10 days to a value 
of about 3000 gauss and remains at that value for about 40 days. Since the growth time is compara- 
tively shorter than the lifetime of the magnetic field we may assume the magnetic field as a function 
of time given by Figure 4 . 

By assuming this magnetic field profile, it is possible to integrate Equation 9. The integration 
yields 


1 . 

2 x 


2 



2 


8x 


- 2 


Ift' - o, x r 1, x = 0 can be assumed as the initial condition, then 

x 2 (l - x 2 ) + y (i “ x ^ + fi (l ” x 2 ) 


( 10 ) 
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MAGNETIC FIELD (gauss) 


where the dot denotes differentiation with respect to t'. Figure 5 shows the relationship of i and x ; 
and x _ l, x x* are the solutions of 


l 



Hi - *-’) 


- 0 . 


3000 

2000 - 
1000 - 

m o 

TIME - 


Figure 4-Profiie used to deduce Equation 10. 




Figure 5-Representation of Equation 19. 


Initially x is an increasing function of t 
because x > o at t “ o. Figure 5 shows 
that x is an oscillatory function of t . A 
complete solution of Equation 9, obtained 
on IBM 7090 computer, is represented in 
Figure 6 for fi - 5, 10, 50. Here, the 
amplitude of oscillation increases with p. 

To obtain a solution of Equation 9, in 
the general case, where the profile of the 
magnetic field is given by Equation 3 we 
must be cautious. If we use the expression 
for IT given in Equation 3 the solution is 
unstable because x > o for all t and for any 
reasonable value of the magnetic field 
strength. For this reason we assume that 
the magnetic field is induced into the area 
at the rate given by Equation 3. The prob- 
lem here is to find a magnetic field profile 
which gives B 0 (t)/x 2 close to the curve 
given in Figure 2 and which gives a 
reasonably good fit for the area variation. 
Assuming then that ^ * 4* * 
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v*) 


l(l - >‘" 2td )(l - ’ I0 ). 


where t is now measured in seconds and d is the number of seconds in one day, we have-instead 
of Equation 9 


d 2 x 

dt 7 



(ii) 


where 


t 


i 



B 2 

g77 Po * 


A solution of Equation 10 for £ " 5, t l 1 is presented in Figure 7. The radius of the cylinder 
reaches a value of 1.6, oscillates about it and finally returns to its original value when the magnetic 
field has disappeared. The period of oscillation is 3 days and the amplitude of oscillation is about 
0.25 r 0 . Figures 8 and 9 present the velocity dx/dt and the temperature ratio T/T 0 as a function 
of t. 
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The assumption^ = 1 corresponds to M/2mip 0 = (86.400) 2 ora = 6 X 10 3 where we have used 
r 0 = 10 9 cm and p = 10 4 . Therefore, a fraction 1/(6 x 10 3 ) of the total mass of the cylinder is 
contained in the current sheath. For much higher values of a, the value of tj becomes large and 
the period of oscillation of the cylindrical sheath may become much smaller than that of Figure 7. 


THE SNOW PLOW MODEL 


In the snow plow model, first discussed by Rosenbluth (Reference 11), the magnetic field is 
homogeneous throughout the interior of the cylinder. The current sheath is in an infinitely thin layer 
of negligible mass. As the sheath expands it collects all the mass with which it comes in contact, 
as a snow plow collects the snow in its path. If r 0 is the initial radius of the cylinder and r the 
radius at any time t, the mass with the sheath is ttNui. ^r 2 - r 0 2 j where n is the average number 
density of ions of mass m i . Newton's equation of motion of the sheath then becomes 


_d_ 

dt 


ttNuIj (r 2 


r 


2 \ Al 

0 ) dt 



Using the transformations 


r 


r 


o 


ft = 



Nm , r o 2 2 
2 Po 


t’ 
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we obtain the equation 



(12) 


where Equations 7 and 8 were used to express B i in terms of B 0 and p L in terms of p 0 . With the 
magnetic field profile of Figure 4, this equation was solved on an IBM 7090 and its solution is shown 
graphically in Figure 10. The initial conditions satisfied by this equation are t ' ■ 0 , x ~ 1, x 
The solution is shown in Figure 10 for p = 5, 10, and 50. 

A peculiar feature of this differential equation is that while in the expanding phase the sheath collects 
all mass with which it comes in contact; however, in the contracting phase it is losing its mass- 
not exactly in the same way as a snow plow because the snow plow loses its mass instantaneously 
when the direction of the velocity is reversed. We can therefore only follow the calculation of the snow 
plow model until the velocity reverses in sign. For that reason this model will not be discussed in 
detail. 



t' 

Figure 1 0— The solution of Equation 12 for 8 = 5, 10, 50. 


THE EFFECT OF DENSITY AND PRESSURE GRADIENTS ON 
THE DYNAMICS OF SUNSPOTS 

The pressure and density both rise in the lower levels of the photosphere. Figures 11 and 12 
show the graph of the variation of the density P and the pressure p from the surface of the photo- 
sphere to the depth of a few hundred kilometers. In the solar photosphere, the pressure rises more 
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Figure 11— Variation of p versus height in the solar 
photosphere (Taken from Reference 10, p. 127, 
Table II). 


300 



>500 l 1 1 1 L J 

icr 10 io ~ 9 io -8 io ~ 7 io ~ 6 io -5 

DENSITY (gm/cm 3 ) 

Figure 12-Variation of p versus height in the solar 
photosphere (Taken from Reference 10, p. 127, 
Table II). 


rapidly with depth than the density. For example at 
a depth of about 400 km the pressure is about 20 
times that at the surface of the photosphere and the 
density is about 10 times that at the surface. If the 
photosphere is assumed to be stratified, the lower 
layers have less amplitude of oscillation than those 
at the surface. Also if the lower layers oscillate 
independently the frequency of oscillation can be 
different at the lower levels. Replacing p 0 with 20p 0 
in Equations 5 and 8, we obtain in place of Equa- 
tion 9: 


I d 2 x 
x dt' 2 



(13) 


The solution of this equation is shown in Figure 13. 



t' 


The frequency of oscillation of x in Figure 13 
is about .75, and in Figure 6 it is about 3.2. Thus 


Figure 1 3— Representation of Equation 13 
for /3 = 5, 10, 50. 
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the two periods of oscillation are 


3-2 y 2wap 0 J ’ 

(■” T af>\ 

1 2 - 0. 3 I 2Trap 0 I 

Therefore the ratio. 75 3.2 = [p'/pV* is greater than or less than one depending on whether p' < 18.2/ 

In order to draw any conclusions about the frequency of oscillation we must know the variation of 
density with depth. From the figures given by Minneart (1953), in Figure 12, we find that p ~ 10 p . 
We therefore conclude that the period of oscillation at lower levels in the solar photosphere is approxi 
mately the same as at the surface. This subject will be taken up in a later paper. 


THE BRIGHT RING AROUND SUNSPOTS 

During the initial growth of the 55-day spot considered the velocity of the current sheath has 
a value of about 1-10 km/sec. Charged particles with their thermal velocities suffer reflection 
from these initially outward moving sheath and gain energy. These particles will then move a few 
mean free paths before loosing the energy gained from reflection. Therefore additional light in the 
solar photosphere or chromosphere, will be generated and a bright ring around a spot, with a width 
of the order of a mean free path will be created. 

The self collision time of particles of mass m and charge e is given by 

1 1 

m 2 ( 3k T) 2 

* c 8 x 0,7147me 4 z 4 log A 


where 


A 


3 / k 3 T 3 \ 2 

2e 2 z 2 \ 7711 / 


Thus the mean free path is 



(3kT) 2 

8 x 0.71477ne 4 z 4 log A 


1.8 x 10 5 T 2 
N log A 


and is independent of the mass of the particle. In Table 1 this quantity is given for N 10 6 to 
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Table 1 


Mean Free Path of a Charged Particle for the Temperatures and Number Densities Shown. 


Temperature 

(°K) 

Mean Free Path of Charged Particles 

N = 108 particles 
cm 3 

N _ 10 g particles 
cm 3 

|*| = ]n 12 particles 
cm 3 

|*| = i ni 5 particles 
cm 3 

10 3 

1.4* 10 4 

19 



! 0 4 

1.1 * 10 5 

1.4 * 10 3 

1.9 


10 5 

.91 * 10 8 

1.1 * 10 5 

1.4 * I0 2 

.19 

10 6 

VI 

'O 

X 

o 

o 

.93 x 10 7 

1.1 * 10 4 

14.5 


10 15 particles/cm 3 and T = 10 3 to 10 6 °K. From this table it appears that the upper chromosphere is 
the region where the mean free path is measurable by optical methods. Near the photosphere the mean 
free path is too short and the thickness of the bright ring will probably be too small to be measured. 

A radially outward moving current sheath moving with a velocity of the order of one km/sec pro- 
duces a shock wave the thickness of which is of the order of 

10 -f- - - 3 ...? 

P (4 X 10 1S X 25 X 1CT 20 ) 2 

- JL 

10 cm 


(Reference 13) which is also too small. Table 2 shows the 
variation of c/&> with density. 

It is only in the upper chromosphere that the thickness 
of the shock is sufficient to be observable. Therefore, ac- 
cording to the present analysis, the bright ring should be 
observed to be expanding. 
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Table 2 


Variation of c/tj p with Number Density. 


N 

/ParticlesN 

c/ % 

V cm 3 / 

10 3 

1.7 x 10 4 

10 6 

540 

10 9 

17 

10’ 2 

.54 

10 15 

7 

o 

X 

IV 
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